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MEASUREMENTS OF RADIATIVE HEATING ON SHARP CONES 

By Dale L. Compton and David M. Cooper 
Research Scientists, NASA, Ames Research Center, Moffett Field, California 


Measurements of total radiative heat load have 
been made on 30 ° half -angle, 1-cm base diameter, 
circular cones in free flight. Models were gun- 
launched into air in the prototype of the Ames 
Hypervelocity Free -Flight Facility with a nominal 
velocity of 7 km/sec. Free-stream pressure was the 
primary test variable and was varied from 2 to 
200 mm Hg. Secondary test variables were angle of 
attack and model material. Tests were performed at 
seven pressures on ablating polyethylene models, at 
four pressures on ablating polyformaldehyde models, 
and at two pressures on nonablating metal models. 

The radiation from the shock layer was measured 
with several broadband radiometers consisting of 
S-5 spectral response photomultipliers. In addi- 
tion, the shock-layer radiation was spectrally 
resolved into 12 components in the 0.22 - 1.15 
micron range with photomultiplier, optical filter 
radiometers . 

Predictions of the total radiation were made 
for both equilibrium and nonequilibrium conditions 
in the shock layer. The equilibrium radiative 
intensities were taken from the work of Kivel and 
Bailey. The nonequilibrium radiative intensities 
were taken from the measurements of nonequilibrium 
radiation of Page and Arnold combined with the 
radiative relaxation distance measurements of Teare, 
Georgiev, and Allen. Radiative intensities were 
computed along streamlines in the conical shock 
layer, and integrations of the local radiative 
intensities were then performed to obtain the total 
radiative heat load. 

Results from the tests show that the radiant 
heat load to cones is much smaller than that to 
blunt bodies at comparable free-stream conditions, 
thus supporting the conclusions of Allen, Seiff, 
and Winovich. 

Radiant power measured at the same free-stream 
conditions by the broadband radiometers was found 
to depend on model material. The radiation from 
the shock layer of the ablating polyethylene models 
was always higher than that from nonablating metal 
and ablating poly formaldehyde models. The addi- • 
tional radiation on the polyethylene models is the 
emission from heated products of ablation in the 
boundary layer, as shown by comparisons of spectra 
for polyethylene with those for metal. Spectra 
from the poly formaldehyde models showed that 
ablation -products radiation was not present in the 
spectral range of the broadband detectors. In addi- 
tion, good agreement was found between the broad- 
band measurements on metal models and those on 
polyformaldehyde models, indicating that for the 
polyformaldehyde models the radiant emission in the 
S-5 spectral response range was primarily from air. 

Comparison of the data from the metal and 
polyformaldehyde models with the predictions of 
radiation for air shows reasonable agreement. At 
the highest test densities , the measured radiation 
is approaching the equilibrium prediction, while 
at lower densities the measurements are substan- 
tially above the equilibrium prediction and agree 


reasonably well with the prediction for combined 
equilibrium and nonequilibrium radiation. 


Introduction 


For entries into the earth's atmosphere at 
velocities well in excess of -escape velocity, Allen, 
Seiff, and Winovich (ref. 1) have shown that the use 
of sharp entry bodies, such as pointed cones, will 
reduce the total aerodynamic heating. This is true 
since the radiative heating which dominates heat 
transfer to blunt bodies at these high velocities 
is reduced to a manageable level on cones by virtue 
of the relatively low gas temperatures behind the 
swept shock wave, and the reduction in radiative 
heating more than compensates for the associated 
increase in convective heating. To help verify 
these conclusions, experimental information is 
required on radiative heat transfer to cones . In 
addition, methods are required for the computation 
of radiative heat load which will account for both 
equilibrium and nonequilibrium radiation in conical 
flow fields. One purpose of the present paper is 
to present the results of free -flight range meas- 
urements of radiative heat load to conical models. 

A second purpose is to develop a method of predict- 
ing the radiative heat load and to compare the 
measurements of radiation with the predictions. 
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Symbols 

area, cm 2 

radiative intensity, watts /cm 3 

equilibrium radiative intensity, 
watts / cm 3 

peak nonequilibrium radiative intensity, 
watts /cm 3 

exponential radiative intensity decay 
factor, cm -1 

cone slant length, cm 

free-stream pressure, 

equilibrium radiating 

free-stream velocity, 

radiant power, watts 

equilibrium radiant power, watts 

Cartesian coordinates 

coordinates defined in sketch (c) 

distance defined in Eq. (B5) 

angle of attack, deg 

angle between the cone surface and 
equilibrium conical shock wave, deg 


mm Hg 

volume, cm 3 
km/sec 


G'Aq /& 7 
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A shock standoff distance in conical layer 

6 p -radiative excitation distance, cm 

6 r distance required for radiation to decay to 
1 . 1 E E q , cm 

0 C cone half -angle, deg 

o y viewing angle defined in sketch (c) 

cp, r\ angles defined in sketch (b ) 

X angle -of -attack correction factor 


Superscripts 


* distances and exponential radiative decay 

factors along a streamline behind a normal 
shock wave in one -dimens ional flow . 


distances along a streamline in a conical 
shock layer 



Models were launched from a light -gas gun into 
the prototype of the Ames Hypervelocity Free -Flight 
Facility, which for this investigation, was 
operated without its countercurrent air stream as 
an enclosed ballistic range. A schematic view of 
the facility is shown in Fig. 1. The model launcher 
is a caliber 0 - 50 , deformable -piston, light -gas gun. 
Models fly through an instrumented test section 
equipped with 11 spark shadowgraph stations spaced 
at 4 -foot intervals along its length. An enlarged 
view of part of the test pection is shown in the 
inset to Fig. 1. Time -of -flight intervals between 
stations are recorded with electronic chronographs. 
The time -distance history and the angle -of -attack 
history of the model are extracted from the shadow- 
graph pictures and chronograph records. Free- 
stream properties are determined from measurements 
of the pressure and temperature in the tunnel just 
prior to launching the model. 

Radiometers 

Radiation from the shock layer of the models 
was measured with six identical broadband radiom- 
eters, depicted schematically in Figs. 1 and 2. 

Each radiometer consisted of an S-5 spectral- 
response photomultiplier (type RCA 1P28). These 
tubes were mounted in housings on the test section 
walls with slits which limited their fields of 
view, and were placed in pairs on orthogonal view- 
ing lines. Thus, each pair of instruments made 
simultaneous measurements on the model from differ- 
ent viewing directions. The spectral range covered 
by the broadband radiometers was from 0.2 to 0.6 
micron between the 10 percent of maximum response 
points . 

Radiation from the shock layer was resolved 
into 1 2 spectral components by another set of 
radiometers also shown schematically in Fig. 1. 

The spectral range covered was from 0*22 to 1.15 
microns. Each radiometer consisted of a photomul- 
tiplier and a narrow -pass optical filter. The 12 
channels of this int rumen t were mounted together 


and viewed the model from nearly the same direction. 
Similar radiometers are described in more detail 
in reference 2 which also describes the methods 
used to calibrate the radiometers for their abso- 
lute radiant sensitivity. Output signals from the 
radiometers were displayed on oscilloscopes and 
photographed. 

Models and Sabots 

The models used for these tests were 30° half- 
angle, right -circular cones with a base diameter of 
1.02 cm. The models were machined from several 
different materials . Plastic models , which ablated 
during flight, were machined from polyethylene and 
from polyf ormaldehyde . Metal models, which did not 
begin to ablate before passing the radiometers, 
were constructed from aluminum with copper tips . 

A sabot machined from polycarbonate plastic was 
used to fit each model to the launch tube of the 
gun. A sketch of a model and sabot is shown in 
Fig. 3- 

Test Conditions 

•The nominal test velocity was 7*0 km/sec. 
However, the actual test velocities varied from 
5*7 to 7.3 km/sec. Polyethylene models were tested 
at free -stream pressures of 2, 5> 10, 20, 50 , 100, 
and 200 mm Hg. Polyf ormaldehyde models were tested 
at pressures of 20, 50, 100, and 200 mm Hg, and 
metal models were tested at pressures of 10 and 
20 mm Hg. Free -stream temperature was room tem- 
perature for all tests. 

Reduction of Broadband Radiometer Data 

The final output desired from the broadband 
radiometer data was the number of watts emitted by 
the conical shock layer of a model flying at zero 
angle of attack and at a velocity of 7 km/sec. To 
obtain this information, the radiometer data were 
reduced in several steps. 

1. The output of the radiometer was multi- 
plied by its calibration constant which was 
adjusted by the inverse square law to correspond 
to the distance from the radiometer to the model. 

2. A correction was then applied to account 
for the difference between the spectral energy 
distribution of the air radiation in the shock 
layer of the model and the spectral sensitivity of 
the photomultiplier. The radiometer responded to 
approximately 80 percent of the air radiation. 

3 . A geometric correction was applied to 
account for the radiation in the part of the shock 
layer hidden from the photomultiplier. At zero 
angle of attack, the radiometer can see approxi- 
mately 70 percent of the shock-layer volume. This 
shadowing effect is relatively insensitive to angle 
of attack. 

4. A correction was applied to account for 
angle of attack of the model. On cones, shock- 
layer temperature and, thus, radiation are sensi- 
tive to angle of attack. To make the correction 
the temperature and density in the conical shock 
layer at angle of attack were assumed to be those 
for the appropriate tangent cone. The predictions 
of reference 3 were then used to obtain the equi- 
librium radiative intensity distributions around 
the cone. These distributions were then integrated 
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over the portion of the model visible to the photo- 
multiplier and divided by the zero angle -of -attack 
value to obtain a correction factor. Corresponding 
confutations were made for nonequilibrium radiation 
in the shock layer. The methods used are presented 
in more detail in appendix A. The success of the 
corrections can be inferred from their effect on 
the data and is discussed in the section on results 
and discussion. 

5. A correction was made to shift the data to 
a standard velocity of 7 km/sec. Each data point 
was multiplied by the ratio where 

is the equilibrium radiant intensity at the test 
condition and Eg^ is the equilibrium intensity 
at the same density but for a free -stream velocity 
of 7 km/sec. The equilibrium intensities were 
again taken from the predictions of reference 3. 

In summary, it was necessary to correct the 
raw data for viewing distance, spectral sensitivity 
of the photomultiplier, geometry of the shock layer, 
angle of attack, and slight differences in test 
velocity to obtain the final data. 

Reduction of Spectral Data 

Each output from the 12 -channel radiometer was 
reduced to watts emitted per micron of wavelength 
by the method described in reference 2. No attempt 
was made to correct these data for angle of attack, 
or to compensate for test velocities different from 
7 km/sec. The principal use for the spectral infor- 
mation was to discover which of the broadband 
radiometer data were being influenced by the pres- 
ence of radiation from heated products of ablation 
in the boundary layer. The influence of this extra- 
neous source of radiation on the data will be con- 
sidered in the section on results and discussion. 


Radiation Predictions 


seen that the region behind the corner contributes 
significantly (approximately 20 percent) to the 
total volume. Equilibrium radiative intensity 
predictions were taken from reference 3. 

Nonequilibrium Radiation 

At low free -stream densities substantial por- 
tions of the shock layer were out of equilibrium, 
and the radiation calculations were extended to 
include nonequilibrium effects. The experiments 
reported in reference 4 show that the nonequilib- 
rium radiative intensity behind a normal shock 
wave has the form shown in Fig. 4(a) - a rapid rise 
to a peak followed by a relatively slow decay to 
the equilibrium level. Since the chemical proper- 
ties and population of excited states in the non- 
equilibrium zone are primarily due to binary 
collisions, the peak of the radiation is approxi- 
mately proportional to the first power of density, 
and the length of the nonequilibrium zone is 
approximately proportional to the inverse first 
power of density. These facts lead to the conclu- 
sion that, as long as the peak intensity is much 
greater than the equilibrium intensity, the heat- 
transfer rate to a body placed at the termination 
of the nonequilibrium zone is nearly independent 
of density for a given free -stream velocity. This 
heat -transfer rate has been measured for some free- 
stream velocities and is reported in references 4, 

5, and 6. In addition, the excitation and relax- 
ation distances (distances for the radiation to 
rise to its peak and to decay to 1.1 times the 
equilibrium value, respectively) have been reported 
in reference 4. For the purposes of the present 
computations, the radiation distribution shown in 
Fig. 4(a) has been approximated by the distribution 
shown in Fig. 4(b) - a linear rise to the peak 
value, followed by an exponential decay to the equi- 
librium level. This approximation gives the follow- 
ing equations for the radiative intensity behind a 
normal shock wave: 


Predictions of the radiation to be expected 
from the conical shock layer for the test conditions 
were made by the following methods. 

Equilibrium Radiation 

To predict equilibrium radiation, the temper- 
ature, density, and volume of the conical shock 
layer must be known. State properties were assumed 
to be constant in the shock layer. To locate the 
equilibrium shock -wave angle, the continuity equa- 
tion for the mass of air swept out by the cone was 
solved simultaneously with the conservation equa- 
tions for the inclined shock wave. This procedure 
is described in detail in appendix B of reference 1. 
The radiating volume was assumed to terminate with 
the first Mach line from the expansion fan at the 
base of the cone. (Radiation from the equilibrium 
expansion fan was also computed and found to be 
small. ) A section through the radiating volume is 
shown in the following sketch, from which it can be . 
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With the additional assumption that the peak radi- 
ation occurs at one tenth of the total distance to 
the point at which the radiation has decayed to 
1.1 E E q, that is, 5p* = 0.1 5 r * (ref. 4), integrat- 
ing Eqs . (l) will give the heat -transfer rate from 
the nonequilibrium layer to a plane at the end of 
the relaxation zone. 
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Performing the integration gives 
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I . t -as v , . .-*"»*) (3) 


In Eq. (3) the unknowns are k*, the exponent 
controlling the decay rate, and Ep, the peak value 
of the radiative intensity. 

A second equation relating Ep and k* is 
obtained when Eqs . (l) are written for x* = 5j.* 

-9k*5p* 

!.! Ejjq = (Ep - EgQ)e + Ej-q (4) 


Solving Eqs. (3) and (4) simultaneously gives 
values for Ep and k* as functions of density p^. 
These equations are valid .when Ep is much greater 
than EgQ, that is, at low densities. For higher 
densities, where this inequality no longer holds, 
Epp was assumed constant in accord with binary 
scaling, and the value of k* was again computed 
from Eq. (4). The results of the computations are 
shown in table I and apply to the streamlines 
behind a normal shock wave at a velocity of 
4.0 km/sec. This velocity was chosen as the veloc- 
ity normal to the conical swept shock wave for a 
30° half -angle cone at a free -stream velocity of 
7 km/sec. Relaxation times were computed from the 
relaxation distances by assuming that the velocity 
behind the normal shock wave changed exponentially 
from the full frozen value immediately behind the 
shock wave to 1.01 times the equilibrium level at 
the time when the radiant level had reached 1.1 
times the equilibrium level. 


To apply these results to cones, it was 
assumed that the relaxation times were the same in 
the conical layer as behind the normal shock wave 
for corresponding free -stream densities and shock 
strengths. Since the air velocity in the conical 
layer is much higher (approximately a factor of 10) 
than behind the normal shock wave, the air flows 
approximately 10 times farther along the cone before 
the radiation returns to its equilibrium level. 
Therefore, the excitation and relaxation distances 
and the exponential radiative intensity decay fac- 
tors were adjusted to account for this high tangen- 
tial velocity. The adjusted values for the conical 
shock layer are also given in table I. It should 
be noted that the radiant heat -transfer rate from 
the nonequilibrium zone to the surface of the cone 
is the same as the heat -transfer rate to a body 
placed at the termination of the nonequilibrium 
zone behind a normal shock wave in one -dimensional 
flow. 


With the radiant intensity distribution along 
streamlines in the conical layer known, the total 
radiant output can be obtained by integrating the 
local radiant intensity over the volume of the 
shock layer. Two cases arise: 

1. The shock layer is near equilibrium. For 
this case, the shock -wave shape was assumed to 
differ negligibly from a cone with the equilibrium 
shock angle. The integration then gave a closed- 
form solution for the heat -transfer rate to the 
conical surface and for the total radiant output 
of the shock layer. 

2. The shock layer is substantially out of 
equilibrium. For this case, the shock wave is no 
longer straight, and the shock -layer volume is 
larger than for equilibrium. The location of the 
shock wave was confuted by assuming the density 


distribution along streamlines to be exponential in 
the same manner as the velocity distribution 
described earlier. The integration of the radiant 
intensity was then performed numerically. Details 
of these procedures are described in appendix £. 

A special problem arises when the density is 
low enough that the nonequilibrium radiative 
profile is severely truncated; that is, when the 
radiation from the streamline entering at the tip 
of the cone does not reach its peak value before 
the stre aml ine reaches the base of the cone. When 
this occurs, the calculated radiative power becomes 
quite sensitive to the manner in which the radia- 
tion is assumed to rise to Ep. As can be seen 
from Fig. 4(a), the assumption of a straight line 
rise to Ep may not be the best choice. In order 
to pick a rise with zero slope at 5p, a quarter 
cycle of a sine wave was used. However, compari- 
son with the actual rise histories from reference 4 
show that the sine wave somewhat overpredicts the 
radiant intensity and thus will overpredict -the 
total radiant power. Predictions for both a 
straight line rise and a quarter sine wave rise 
are shown in the final results. 

The procedure described above is open to crit- 
icism on several points. First, it does not take 
into account any viscous effects; calculations by 
the methods described in reference 7 show that for 
the lowest test densities, the displacement thick- 
ness of the boundary layer is an appreciable frac- 
tion of the shock -layer thickness even for the 
nonablating models. The thick boundary layer will 
tend to push the shock wave away from the cone, 
effectively increasing the shock strength. and thus 
increasing the radiant intensity. Furthermore, the 
high -temperature air in the boundary layer will 
contribute to the radiation. Second, there are no 
experiments to verify either the equilibrium radi- 
ative intensities predicted for these temperatures 
and densities or the nonequilibrium radiative heat- 
ing rates, since the conditions for the present 
tests correspond to normal shock velocities below 
those previously studied. In fact, the nonequilib- 
rium heating -rate data from reference ? must be 
extrapolated from a velocity of 6 to 4 km/sec and 
this extrapolation is certainly open to question. 
Third, no attempt was made to treat the expansion 
wave at the base of the cone under nonequilibrium 
conditions. These difficulties, along with the 
other approx ima tions already noted, make the method 
approximate at best. 

In spite of the disadvantages, the correlation 
of the predictions w-ith the data is reasonable, as 
will be shown later. The method has advantages as 
well. It is relatively simple to use and does not 
require detailed computations of the local thermo- 
dynamic properties in the conical shock layer. It 
can be used to predict local radiant heat -transfer 
rates to the conical surface as well as the total 
radiant power in the conical shock layer. 


Results and Discussion 


The test results demonstrate experimentally 
the validity of the thesis of reference 1 - that 
radiative heat load is vastly smaller on cones than 
on comparable blunt bodies. In fact, the radiative 
heat lord measured in the present tests on 30° 
half -angle cones is 2 to 3 orders of magnitude less 
than the radiative heat load reported in reference? 
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on "blunt "bodies (when adjustments of the data 
from reference 5 are made to comparable test con- 
ditions and "body size). This reduction in radia- 
tive heat load results directly from reduced shock- 
layer temperatures, which, in turn, result from 
substitution of oblique conical shock waves for 
strong detached shock waves. 

The broadband radiometer data from the present 
tests are shown in Fig. 5, where the radiant power 
emitted from the shock layer in watts (corrected 
for the effects enumerated in the section on data 
reduction) is plotted as a function of the calcu- 
lated equilibrium shock-layer density in amagats. 
Each point plotted on Fig. 5 represents an average 
of the data points obtained for that density and 
model material. Data are shown for the three model 
materials: polyethylene, metal, and polyformalde- 

hyde. It can be seen from the figure that the 
radiation from the ablating polyethylene models is 
always greater - by factors of 1.7 to 4 - than that 
from the nonablating metal models and the ablating 
poly formaldehyde models. This is due to radiation 
from heated products of ablation in the boundary 
layer on the polyethylene models. Spectra from a 
polyethylene model and a metal model at identical 
free -stream conditions are shown in Fig. 6. It 
can be seen that the spectrum from the ablating 
polyethylene shows radiation in the visible and 
infrared which exceeds that in the spectrum from 
the nonablating metal model. The spectral range 
covered by the broadband detectors is also shown, 
and it is evident that it includes part of the 
radiation from ablation products. Thus the results 
shown in Fig. 5 for the polyethylene models are 
influenced by this extraneous source of radiation. 
Ablation -products radiation has been studied exper- 
imentally for blunt bodies in reference 2, where a 
method was developed for its correlation, and for 
cones in reference 7; in which the correlation 
parameters from the blunt -body tests were found to 
predict successfully the measured ablation -products 
radiation on cones . Reference 7 also points out 
that radiation from the hot surface of the model is 
not sufficient to account for the observed radia- 
tion, as is demonstrated on Fig. 6 where the maxi- 
mum possible surface contribution is shown 
multiplied by 100 for an assumed wall temperature 
of 1000° K. It was not possible to obtain data on 
nonablating metal models at higher densities with 
the present equipment, since the extremely high 
heating rates at the sharp tip of the model would 
cause ablation to begin before the model reached 
the radiometer stations. Tests were therefore 
conducted with polyformaldehyde models for which 
the spectra (unpublished) from blunt -body tests 
indicated that radiation from ablation products 
would negligibly influence the broadband measure- 
ments. The resulting spectra on cones also showed 
that the radiation from ablation products in the 
range of the broadband detectors was small. The 
good agreement shown on Fig. 5 between the results 
obtained with the ablating polyformaldehyde models 
and the nonablating metal models at 0.25 amagat 
further supports this conclusion. Therefore the 
radiant power shown on Fig. 5 for the polyformalde- 
hyde models is primarily from air. 

Fig; 5 also shows the predictions for air 
radiation in the conical shock layer made by the 
methods described earlier. Predictions are shown 


for both equilibrium and combined equilibrium- 
nonequilibrium radiation. The agreement between 
the measurements made on metal and polyformaldehyde 
models and the predictions which include both equi- 
librium and none qui lib rium effects is reasonably 
good. At the highest densities the measurements 
are approaching the equilibrium predictions. The 
measurements made on polyethylene models are sub- 
stantially above the prediction at all densities, 
primarily because of the influence of ablation - 
products radiation discussed earlier. However, it 
should be noted that at the lower test densities, 
viscous effects may also be important, since at 
these low densities the boundary -layer thickness 
is becoming a substantial fraction of the shock- 
layer thickness. 

All of the data in Fig. 5 are presented as if 
the cones were at zero angle of attack. Since this 
was never actually the case, corrections were 
derived to account for the resulting change in 
radiant power. The method for deriving these cor- 
rections is presented in the section on data reduc- 
tion and in appendix A. If these corrections are 
successful, they should have two effects on the raw 
data. The first is to give the same deduced radia- 
tive power for each pair of radiometers placed at 
the same measuring station, but viewing the model 
from different directions; and the second is to 
make the deduced radiative power appear independ- 
ent of angle of attack. 

Typical results are shown in Fig. 7 for the 
polyethylene models at a free -stream density of 
100 mm Hg, where the radiant power is plotted as 
a function of angle of attack. (These data have 
already been corrected for slight velocity differ- 
ences.) The circles show the data before the 
angle -of -attack correction was applied and the 
squares show the data after correction. Pairs of 
data points at the same a are from simultaneous 
measurements. It can be seen that, for the most 
part, the corrected data have the desired behavior, 
indicating that the corrections are in fact reason- 
ably successful in removing the angle -of -attack 
effects . 

It will be noted in Fig. 7 that the uncor- 
rected data do not continue to Increase with 
increasing angle of attack. This is because the 
radiometer does not always view the model from the 
same direction (see sketch (c), appendix A). Thus 
the measured radiation depends on both angle of 
attack and radiometer viewing angle. The angle - 
of -attack correction was derived to account for 
both of these effects. 

After the corrections had been applied to the 
broadband data, considerable scatter existed for 
which causes could not be found. Root -mean -square 
deviations from the average values of radiant power 
are given in table II. No root -mean -square devia- 
tion is given for the polyethylene model data at 
a free -stream pressure of 2 mm Hg because of insuf- 
ficient ' data.’ 


Conclusions 

From the results of these studies of radiant 
emission from the shock layer of conical models 
it can be concluded that: 
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1. The measurements of radiation show that 
the substantial reduction in radiative heat load 
predicted for sharp cones over blunt bodies is, in 
fact; realized. 

2. A method has been developed for predicting 
in some detail the radiative heat load and heat- 
trains fer rate to cones from air in the shock layer. 
The method accounts for both equilibrium and non- 
equilibrium radiation. 

3. Measurements of radiation on nonablating 
metal models emd ablating polyformaldehyde models 
compared with the predictions for air radiation 
show reasonable agreement. 

4. Measurements of radiation on ablating 
polyethylene models show that the radiation from 
heated products of ablation in the boundary layer 
provides a net increase in heat load, as shown by 
comparison of data from the ablating polyethylene 
models with data from the nonablating metal models. 

5. The effects of angle of attack on the data 
have been successfully correlated. 


free -stream densities for which the shock layer was 
near equilibrium, this ratio was nearly independent 
of free -stream density. The ratio Veq/(Veq ) (1 ^ 0 
was also computed. This ratio is the relative 
increase or decrease in local shock-layer volume as 
a function of (p and a and was also found to be 
nearly independent of free -stream density. 

These two ratios were combined to give 
e EQ v Eq/ ( e EQ^EQ ) a - Q * This ratio represents the 
local value of the radiation in watts as a function 
of cp and a, divided by the value at a - 0. To 
form a correction factor it was integrated over the 
portion of the conical shock layer visible to the 
photomultiplier. For this purpose it was assumed 
that the photomultiplier viewed 64 percent of the 
volume before the shoulder of the model (this per- 
centage was determined graphically) and all of the 
volume between the shoulder and the first wave of 
the expansion fan. Computations of the total radi- 
ating volume showed that 80 percent of the volume 
was in front of the shoulder and 20 percent was 
between the shoulder and the expansion fan. 

The angle -of -attack correction is therefore 
given by 


Appendix A 

Angle -of -Attack Correct ions 

The following procedure was used to form a 
correction factor to account for the effect on the 
radiometer data of angle of attack of the model. 

The coordinate system is shown in sketch (b) 
where cp is the angle around the base of the cone, 
and cp = 0 is at the ray on the cone surface which 
experiences the maximum inclination toward the free- 
stream .velocity vector. 



State properties in the shock layer, for a • 
given value of cp, were assumed to be constant 
along lines perpendicular to the conical surface 
and were calculated by the mass -flow method, 
described in the text, by assuming that they were 
the properties on the appropriate tangent cone. 
Thus, the distributions of temperature, density, 
and shock -wave angle were assumed to be functions 
only of the angle r\ between rays on the surface 
of the cone and the velocity vector. 


X (c v ,a) 
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where the first term on the right accounts for the 
volume before the shoulder and was integrated over 
the portion of the cone visible to the photo- 
multiplier. The second term accounts for the 
volume between the shoulder and the expansion fan 
and was integrated completely around the cone. The 
viewing angle, a v , which defines the direction from 
which the photomultiplier views the model, is shown 
in sketch (c). 


RADIOMETER 


RADIOMETER 


PROJECTION / 

OF 1 ON X J, | y 
X 2 , X 3 PLANE - 

MODEL <t_* 

VECTOR 


/ 


RADIOMETERS ARE 
IN X 2 , X 3 PLANE 


Sketch (c) 


The relation between cp, a, and q is 

cos q = cos a cos 0 C - cos cp sin a sin 0c (A1 

With temperature, density, and shock -wave angle 
known as functions of <p and a, the ratio 

was as a function of cp 

and a by use of the equilibrium predictions of 
reference 3* It was found that, for the range of 


Figure 8(a) shows the values of X plotted 
as a function of viewing angle for various angles 
of attack. The corrections shown are for equilib- 
rium and, when applied to the raw data, it was 
apparent that, at the lower test densities, the 
corrections were too large. Since the radiation 
at these lower densities was known to be primarily 
nonequilibrium, the corrections were recalculated 
with the nonequilibrium radiation measurements from 
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0 < x* < &p 


reference 5 substituted for the equilibrium values. 
The computed corrections are shown in Fig. 8(b) to 
be considerably smaller than the corrections for 
equilibrium. This is because the dependence of 
radiation on velocity normal to the shock wave is 
much less for nonequilibrium than for equilibrium. 
When the corrections for nonequilibrium were 
applied to the test data at low densities, the 
results were adequate. Equilibrium corrections 
were used for free- stream pressures of 100 and 
200 mm Hg, and nonequilibrium corrections were used 
for pressures of 50 11211 Hg an ^- below. 

It is evident that a straightforward extension 
of this method could be used to obtain estimates of 
the local heat -transfer rates and total radiant 
heat load on cones flying at angle of attack. 

A better approach would be to treat each data 
point individually, by dividing the shock layer 
into a number of streamtubes and computing the radi- 
ation from each streamtube separately, taking into 
account the distribution of radiation along each 
streamtube. However, this approach appeared to be 
prohibitively time consuming and, in view of the 
reasonable results obtained with the simpler meth- 
ods, unnecessary. 


Appendix B 

Integration of Radiation in the Conical Shock layer 

The text describes the models of radiation used 
to obtain the radiation distribution from air along 
streamlines in the conical shock layer. In order 
to obtain estimates of the total radiant power from 
these distributions, it was necessary to perform a 
volume integral of the local radiant intensity over 
the radiating volume. 


The coordinate system used for this purpose is 
shown in the sketch below 



SHOCK WAVE 


STREAMLINE 


CONE SURFACE 


Sketch (d) 

The shock standoff distance is 

A = x tan 7 (Bl) 

For purposes of this integration, streamlines were 
assumed to be straight lines parallel to the body 
surface. The distance along a streamline measured 
from the shock wave is given by x 1 . The y coor- 
dinate of a streamline is constant and is given by 

' y = (x - x*)tan 7 (B2) 

Combining Eqa. (Bl) and (B2) and rearranging gives 


x 1 


A- .Z 

tan 7 


(B3) 


The distribution of radiation along streamlines is 
given by Eqs. (l) of the text, with the notation 
changed to that for streamlines in the conical 
shock layer, as 


E(x') = Ep^ 

— k(x * - 5n) 

E(x') = (E P - E EQ )e ^ + EgQ x* > 6p 


(B4) 

Making the change of variables in Eqs. (b 4) to the 
cone-apex-oriented coordinate system described by 
Eqs. (Bl), (B2), and (B3) gives for the distribu- 
tion of radiation in the conical shock layer 


A - y 
tan 7 


E(x,y) = E P - y = (A - &p tan 7) < y < 
°P 


E(x,y) = (Ep - EpQ)e 


(&r*) 




+ Eeq y < ypj 
(B5) 


Integrating first with respect to y gives 
the infinite, thin, plane shock- layer idealization 
to the local heat- transfer rate. The integrals are 



0 < x < 6p| 


W 

A 



dy > (B6) 



x > &p J 


where w/A is twice the local heat-transfer rate 
and the area A is a function of x. This factor 
of 2 occurs since the equations in this appendix 
are written for the total radiant emission from 
the shock layer, of which only about 50 percent is 
incident on the surface. Integrating Eqs. (B 6) 
gives 


W = Sp x£ 
A tan 7 2 Sp 


0 < x < &p| 


A tan 




>(B7) 


+ ^ &p 
2 


x > &p 


Integrating once more, over the surface area 
of the cone, will give the total radiant power. 

The surface area element is given by 

dA = 2rt x sin 0 C dx (b8) 

where 0 C is the half apex angle of the cone. 
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Thus the power can he written 



It is convenient to normalize Eqs. (B9) with 
respect to the total radiation that would be 
obtained if the shock layer were fully in equilib- 
rium. For equilibrium, the radiation along stream- 
lines is constant and twice the heat-transfer rate 
is given by 

Tp = e eq 4 (bio) 

The total equilibrium radiative power is given by 
the integration of Eq. (BIO) with respect to x 
to give 

W E q = 2it tan 7 sin 0 C J E^qX 2 dx (Bll) 


Integrating Eqs. (B9) and (Bll) and taking the 
ratio W/W E q gives finally 



- (k&p + l)j | l > &p J 


Equations (B12) give the radiant power in the 
conical shock layer for the volume ahead of the 
cone base. They do not include the power in the 
volume behind the cone base and before the expan- 
sion fan. (See sketch (a) for the definition of 


radiating volume.) To partially compensate for 
this difficulty, the value of VTgQ for the total 
radiating volume was used in applying equation (B12) 
to compute the radiation expected in the present 
tests. This assumes that the ratio W/WpQ is the 
same for both volumes - an assumption which is not 
exact. However, the error resulting from this 
procedure is small since the volume behind the cone 
base conprises only 20 percent of the total volume. 

Two other effects have been considered. The 
first of these occurs since, for nonequilibrium 
flow, the shock wave is curved and lies outside the 
equilibrium wave. Thus the integration procedure 
described above, which is over only the equilibrium 
volume, will underestimate the total radiative 
power. The effect will be largest at the lower 
densities. To account for this effect, the density 
was assumed to rise exponentially (along stream- 
lines) from its frozen value immediately behind the 
shock wave to 0-99 of its equilibrium value at the 
point on the streamline where the radiation has 
decayed to 1.1 E E q. The position of the shock 
wave was located by the mass -flow method described 
in the text for the nonconstant density. The 
equations for the total radiant power were then 
integrated numerically. It was found that the 
effect was small for free -stream pressures greater 
than 20 mm Hg. This effect has been included in 
the predictions presented in Fig. 5. 

The second effect is described in the text. 

A sine wave was substituted for the initial lin- 
ear rise of the radiation to Ep. The results are 
also shown in Fig. 5- 
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TABLE I. - CONSTANTS USED IN RADIATION CALCULATIONS FOR TEST CONDITIONS 



Normal shock 
wave 

Conical shock 
layer 

Poo. 
mm Hg 

Ep, 

watts /cm 3 

e eo> 

watts/cm 3 

k ti 

cm 

8p*> 

cm 

k ’-i 

cm 

Sp 

cm 

5 

1.77 

0.116 

3.02 

0.0182 

0.235 

0.234 

10 

3.20 

.297 

5.58 

.00912 

' -435 

• n -7 

20 

4.55 

.970 

8.80 

.00456 

.684 

.0586 

50 

15-9 

4.50 

19.7 

.00182 

1-54 

.0234 

100 

31-5 

13-5 

31.6 

.000912 

2.46 

.0117 

200 

63.0 

40.5 

41.8 

.000456 

3.25 

.00586 


TABLE II. - ROOT-MEAN-SQUARE DEVIATION OF BROADBAND RADIOMETER DATA 


Model 

material 

Free -stream 
pressure, 
mm Hg 

Average watts 
observed 

: (FiS- 5 ) . 

Number of . 
data points 

Root -mean -s quar e 
deviation, 
watts 

Polyethylene 

200 

4.1 

10 

±0.62 


100 

3.9 

12 

i -59 


50 

2.0 

10 

±.53 


20 

.58 

4 

±.12 


10 

•34 

8 

±.12 


5 

.18 

11 

±.045 

Polyformaldehyde 

200 

2.0 

14 

+.26 


100 

1.0 

12 

±.15 


50 

•71 

10 

±.21 


20 

.31 

2 

±.084 

Metal 

20 

.36 

23 

±.13 


10 

.21 

8 

±.063 
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Pig. 5-- Measurements and predictions of total radiant power for 1.02 cm base 
diameter, 30 ° half -angle cones; * T-0 W««c- 


* 100 mm Hg 
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Pl«. 7- Typical effect of angle -of -attack correction on data; V„'= 7*0 km/sec. 




Pig. 9 -- Angle -of -attack correction curves, (a) Equilibrium radiation, (b) 




